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Abstract. We recently reported the photometric and spectroscopic detection of the primary transit of the 1 1 1 -day-period, 
eccentric extra-solar planet HD 80606b, at Observatoire de Haute-Provence, France. The whole egress of the primary transit 
and a section of its central part were observed, allowing the measurement of the planetary radius, and evidence for a spin-orbit 
misalignment through the observation of the Rossiter-McLaughlin anomaly. The ingress having not been observed for this 
long-duration transit, uncertainties remained in the parameters of the system. We present here a refined, combined analysis of 
our photometric and spectroscopic data, together with further published radial velocities, ground-based photometry, and Spitzer 
photometry around the secondary eclipse, as well as new photometric measurements of HD 80606 acquired at Mount Hopkins, 
Arizona, just before the beginning of the primary transit. Although the transit is not detected in those new data, they provide 
an upper limit for the transit duration, which narrows down the possible behaviour of the Rossiter-McLaughlin anomaly in 
the unobserved part of the transit. We analyse the whole data with a Bayesian approach using a Markov-chain Monte Carlo 
integration on all available information. We find R p = 0.98 ± 0.03 Rj up for the planetary radius, and a total primary transit 
duration of 11.9 ± 1.3 hours from first to fourth contact. Our analysis reinforces the hypothesis of spin-orbit misalignment in 
this system (alignment excluded at > 95 % level), with a positive projected angle between the planetary orbital axis and the 
stellar rotation (median solution A ~ 50°). As HD 80606 is a component of a binary system, the peculiar orbit of its planet could 
result from a Kozai mechanism. 



1. Introduction 

HD 80606 is a solar-type star with a gas giant planetary com- 
panion on a highly eccentric 111-day orbit (Naef et al. 2001). 
With e = 0.93, the planet receives about a thousand times 
more star light at periastron than at apastron, which makes 
it a key system to study the atmospheric and thermal proper- 
ties of hot gas giant planets. By a lucky coincidence (about 
1 percent probability for a randomly oriented orbit), the or- 



* Based on observations made with the 1.20-m and 1.93-m tele- 
scopes at Observatoire de Haute-Provence (CNRS), France, by the 
SOPHIE consortium (program 07A.PNP.CONS), and with a 16-inch 
telescope at Mt. Hopkins, Arizona, USA, by the MEarth team. 



bital plane is aligned with the line-of-sight, so that both the 
secondary eclipse and primary transit were detected. The sec- 
ondary eclipse was measured during a long photometric run 
with the Spitzer space telescope (Laughlin et al. 2009). In 
Moutou et al. (2009, hereafter M09), we presented our detec- 
tion of the primary transit, simultaneously measured in pho- 
tometry and spectroscopy with the 1.2-m and 1.93-m tele- 
scopes at Observatoire de Haute-Provence, France. The spec- 
troscopic transit data seemed to indicate that the orbital plane 
of the planet was not aligned with the stellar rotation axis. But 
since the transit ingress was not observed, a large degree of un- 
certainty remained in this parameter, as well as in the latitude 
of the transit and radius of the host star. 
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Photometric data of the same event are available from two 
other teams and locations. Fossey et al. (2009, F09) presented 
data obtained at the Mill Hill London Observatory, England, 
with two telescopes (a 35-cm Celestron and a 25-cm Meade), 
on a time span that almost matches that of our OHP obser- 
vations the same night: data were obtained during the main 
portion of the flat part of the transit, the whole egress, and a 
few hours after its end. Garcia-Melendo & McCullough (2009, 
G09) presented data obtained with a 60-cm telescope at the 
Esteve Duran Observatory, Spain, on a shorter time span: the 
observations started just before the egress. No detection of the 
transit ingress has been reported at the time of writing. 

In this paper, we present photometric data of HD 80606 
taken from Mt Hopkins, Arizona on the same night with the 
MEarth network (Irwin et al. 2009; Nutzman & Charbonneau 
2008). These data show no flux variation, but they do provide 
a powerful constraint on the system parameters by imposing 
a strict lower limit to the beginning of the transit. We apply 
a Bayesian analysis to the whole data set, together with pre- 
vious radial-velocity monitoring and the Spitzer observations 
near secondary eclipse, to calculate accurate values of the sys- 
tem parameters, including the radii of the host star and planet, 
and the spin-orbit angle. 

2. Observations 

2.1. MEarth photometry 

A single field containing HD 80606 and HD 80607 was mon- 
itored continuously using one telescope of the MEarth obser- 
vatory located at the Fred Lawrence Whipple Observatory on 
Mount Hopkins, Arizona, for the night of 2009 February 13. 
Additional observations were taken on 2009 February 14th and 
15th but these are not used in the present work. MEarth uses a 
non-standard 715 nm long-pass filter, with the response limited 
at the red end by the long-wavelength tail of the CCD quantum 
efficiency curve. 

Observations were started at the end of nautical twilight 
(solar elevation 12° below the horizon), and continued until the 
start of nautical twilight in the morning. The airmass of the 
field varied from 1 .9 at the start of observations, to a minimum 
of 1.1 at meridian transit, which occurred at UT 07:11, and at 
the end of observations was 2.5. A total of 1695 4.8 s exposures 
were obtained at an average cadence of ~ 23 s including over- 
heads (CCD readout time, and re-centering the field after each 
exposure). 

Due to the extremely short exposures necessary to avoid 
saturation, and the small telescope aperture of 0.4 m, errors on 
individual measurements are dominated by atmospheric scin- 
tillation. It was not possible to defocus the telescope due to the 
small (~ 21") separation between HD 80606 and HD 80607. 
We use the formula of Young (1967) to estimate the contribu- 
tion to our observational error bars arising from scintillation. 
As noted by Ryan & Sandler (1998), the typical coherence 
length for this effect is ~ 12", so the standard formula should 
remain a reasonable approximation. 

Data were reduced using the standard MEarth reduction 
pipeline, which is at present largely identical to the Monitor 



Date [HJD] 
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2454875.591749 


8.153709 


0.005435 


2454875.592015 


8.074484 


0.005425 



Table 1. Photometric times series for HD 80606 from MEarth (full 
table available electronically) 



project pipeline described in Irwin et al. (2007). An aperture 
radius of 10 pixels (corresponding to 7'.'6 on-sky) was used to 
extract differential photometry. In order to estimate the con- 
tamination in the aperture for one star from the other given this 
large aperture size, we use measurements in multiple concen- 
tric apertures from single-stars in the same field as our target 
to derive a simple curve-of-growth. From this we estimate that 
< 0.3% of the flux from HD 80607 falls inside the aperture 
centered on HD 80606, and vice versa. This is negligible for 
our purposes, since it would lead to a < 1% underestimation in 
transit depths, which is smaller than the observational error in 
this quantity. 

We used HD 80607 as a comparison star to derive differen- 
tial light curves of HD 80606 from this photometry. Since these 
stars have similar positions on the detector and almost identical 
colors, effects such as color-dependent atmospheric scintilla- 
tion and flat fielding error are minimized by doing this, and we 
find no advantage to attempting to use other stars of comparable 
brightness on the field as additional comparison stars. MEarth 
uses German Equatorial Mounts, so the entire telescope and de- 
tector system must be rotated through 180° relative to the sky 
upon crossing the meridian. Using HD 80607 as a comparison 
star, we see little evidence for flat fielding errors in the data 
for HD 80606, which normally manifest as different base-line 
levels in the light curve for positive and negative hour angle. 
We therefore apply no correction for this effect in the present 
analysis. 

The MEarth data is given in Table 1 . 
2.2. OHP 120-cm photometry 

The photometric observations of HD 80606 and HD 80607 per- 
formed at the 120-cm telescope at OHP during the nights 2009 
February 12 and 13 were presented by M09. The transit was 
detected during the second night, which is the only one of the 
two that we use in the present work. 326 frames were secured 
during the transit night, with 20 to 30-second exposure times. 
The negative slope after egress tends to suggest that a correc- 
tion for airmass variations should be taken into account in de- 
riving the flux. The airmass is 1 .5 at the beginning of the night, 
then reaches 1.0 at the middle of the egress, and increases up 
to 1.66 at the end of the night. Assuming that the slope seen 
out of transit is due to airmass changes, and that the effect of 
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airmass on the photometry is linear, we find a correction of 
(2 - airmass) x 0.004. The impact of the correction is to create 
a small slope at the beginning of the transit sequence, i.e. dur- 
ing the flat section of the transit. The significance of this new 
slope is low. Since this correction dominates the error budget 
on the transit parameters, it is important to take correlated noise 
into account in the analysis. 



2.3. SOPHIE Doppler spectroscopy 

Radial velocities of HD 80606 measured with the SOPHIE 
spectrograph at the 1.93-m telescope of OHP were presented 
by M09. A continuous sequence of 39 measurements was ac- 
quired during the transit night (2009, 13-14 February), as well 
as nine extra measurements between 8th and 18th February, 
out of the transit. The radial velocities were obtained from a 
weighted cross-correlation of the spectra with a G2-type nu- 
merical template. All the spectra present a similar signal-to- 
noise ratio, S/N 47 per pixel at 550 nm. Together with the 
parameters of the cross-correlation function of the spectra (full 
width at half maximum of 7.33 + 0.02 kms -1 , and a contrast 
representing 48.2 ± 0.3 % of the continuum), this corresponds 
to a ~ 2.4 ms 1 photon-noise uncertainty on the radial ve- 
locity measurements. We quadratically added 3.5 ms _1 due to 
telescope guiding errors and 1 ms _1 due to wavelength cali- 
bration, resulting in radial velocities with a typical accuracy of 
~ 4.5 ms _1 . 



2.4. Published Doppler and photometric data 

We complemented our SOPHIE measurements with published 
radial velocities from three other instruments: 74 measure- 
ments from ELODIE at OHP (Naef et al. 2001; Moutou et al. 
2009), 46 HIRES measurements from Keck (Butler et al. 2006), 
and 23 HRS measurements from the Hobby-Eberly Telescope 
(Wittenmyer et al. 2009). The typical accuracies of these 
three extra datasets are 12.3, 5.0, and 8.3 ms~' respectively, 
and their time spans are Nov. 1999-Nov. 2003, Apr. 2001- 
Feb. 2005, and Dec. 2004-Mar. 2007, respectively. None of 
those measurements were obtained during a primary transit. 
These measurements are used to refined the Keplerian orbit of 
the planet. Systematic radial velocity shifts between the dif- 
ferent datasets are unknown, and left as free parameters. The 
residuals around the model orbit show that the ELODIE uncer- 
tainties are underestimated, and we scaled these uncertainties 
upwards by a factor 1.8 in order to obtain coherent normalized 
residuals. 

The egress of the transit of 14th February 2009 was ob- 
served in photometry from three locations covering part of 
the transit center and the transit egress (M09, F09, G09; see 
Section 1). 

HD 80606 was also monitored for nearly 24 hours around 
the time of the secondary eclipse at 8 microns with the Spitzer 
space telescope (Laughlin et al. 2009). These data provide 
strong additional constraints on the system parameters by mea- 
suring the time and duration of the secondary eclipse. 



Radial velocity 




0.4 0.6 
Phase 



Fig. 1. Radial-velocity data and model. Data from Elodie (blue), 
Sophie (red), Keck (magenta) and HET (green). 



We include these data sets in our combined solution (we 
use only the Celestron time series from F09). 



3. Analysis 

Figures 1, 2 and 3 display the photometric and spectroscopic 
data for the radial velocity curve, the spectroscopic transit and 
the photometric transit. The secondary eclipse is plotted in 
Laughlin et al. (2009). Figure 4 shows the configuration of the 
system according to the best-fit solution in Section 4. 

The procedures to infer physical parameters from observa- 
tions of transiting planetary systems have now been firmly es- 
tablished, and numerous descriptions can be found in the recent 
literature (see e.g. contributions to Pont et al. 2009, for reviews 
and examples). The case of HD 80606b requires all the toolbox 
of the trade for several reasons: the extreme eccentricity makes 
the relation between physical parameters and observable quan- 
tities even more non-linear than usual, the incomplete coverage 
of the transit in photometry and spectroscopy requires sound 
Bayesian statistics, most data sets are dominated by correlated 
noise rather than random errors, and different types of informa- 
tion (transit and secondary eclipse photometry, radial velocity 
orbit, stellar evolution models) provide partial constraints of 
comparable importance. For these reasons, we need to use a 
fully Bayesian approach, with physically meaningful priors on 
all parameters. We also need proper accounting for correlated 
noise. 

Although the data is incomplete in some respects (absence 
of detection of the transit ingress), the abundant datasets leave 
little leeway for most parameters. In particular, the combination 
of the secondary eclipse duration, the egress duration and the 
orbital parameters tightly constrain the mass and size of the 
host star and the planet. 
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Fig. 2. Radial-velocity data and models around the phase of transit. The three model curves correspond to spin-orbit angles at the best-fit 
solution (A = 33°, solid line), a good solution with a higher spin-orbit angle (A = 122°, dashed), and the highest-likelihood solution with 
aligned spin and orbit (A = 0, dash-dotted). 
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Fig. 3. Photometric data on 13-14 Feb 2009 transit. Red: M09, black: MEarth, blue: F09, green: G09 



3.1. Method 

We use a Bayesian approach to determine the probability distri- 
bution of the physical parameters given all the available obser- 
vations and prior information such as stellar evolution models. 
We perform the integration with a Markov Chain Monte Carlo 
(MCMC) algorithm, using Metropolis-Hastings sampling. This 
is similar to the procedure described for instance in Collier 
Cameron et al. (2007). 

This technique is particularly sensitive to the assumptions 
used on the distribution of measurement uncertainties. MCMC 



integration samples the merit function used to measure the 
agreement between model and observations according to the 
Bayesian prior probability, and as such it cannot be better than 
the merit function used. If a traditional "observed vs. mea- 
sured" sum-of-squares is chosen, the error underestimation due 
to correlated noise that plague straightforward fitting methods 
will also be observed. To avoid biased results and unrealisti- 
cally narrow posterior probability distributions for the output 
parameters, it is essential to use realistic error estimates for all 
data, including systematic errors. In the photometric and transit 
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spectroscopic time series of HD 80606 the dominant source of 
systematics are fluctuations correlated in time (e.g. airmass and 
seeing effects). Pont et al. (2006) discuss the effects of this "red 
noise" and a possible way to integrate it in the merit function. 

3.2. Correlated noise 

For most datasets pertaining to HD 80606, red (i.e. correlated) 
noise is by far the dominant source of uncertainties: (1) the 
SOPHIE coverage of the spectroscopic transit was acquired 
during a single night. It is known that in such conditions, 
SOPHIE data are sensitive to weather-related systematics at 
the level of a few meters-per-second; (2) the photometric times 
series also cover only the transit egress. When only a partial 
transit is observed with ground-based photometry, the correc- 
tion of slow trends is difficult and dominates the error budget; 
(3) the out-of-transit baseline flux of the Spitzer observations is 
made variable at the 10~ 4 level by instrumental effects and by 
the light reflected on the day side of the planet, a factor which, 
from the point of view of parameter determination, is equiva- 
lent to red noise. 

We model the red noise as described in Pont et al. (2006), 
with a single 07 parameters for each data set, describing the 
amplitude of correlated noise over the relevant timescale. For 
the photometric time series we use the "V(«)" method de- 
scribed in that paper, which estimates cr r from the departure 
of binned versions of the data compared to purely white noise. 
We find 07 = 7 x 10~ 4 for the OHP photometry, 8 x 10~ 4 for the 
MEarth photometry, 3.6 x 10~ 4 for the F09 photometry , 8x 10~ 4 
for the G09 photometry, 4 ms" 1 for the Sophie transit spec- 
troscopy and 2 x 10~ 4 for the Spitzer time series. The first three 
values are coherent with our experience of ground-based transit 
spectroscopy, with 07 between 4 and lOx 10~ 4 being typical for 
single-object rapid cadence time series. The value of the corre- 
lation parameter for the SOPHIE time series is also compatible 
with previous experience (e.g. Hebrard et al. 2008). Finally, 
there is enough out-of-transit data in the Spitzer time series to 
measure 07 precisely. Our red noise analysis thus indicates that 
the F09 Celestron data is subject to lower systematics than the 
M09 and G09 photometric measurements (as visible as well in 
the behaviour of the residuals compared to model lightcurves). 
This implies that the first data set will get about twice as much 
weight as each of the other two in the combined analysis. To 
account for the clearly much lower quality of the photometric 
data at the beginning and end of the night, we increase 07 by a 
factor three for these measurements (JD < 2454876.34). 

3.3. Models 

We use the Mandel & Agol (2002) algorithm to build model 
transit and secondary eclipse lightcurves, with a linear limb- 
darkening law with u = 0.66 (suitable for a solar-type star; 
at this point the accuracy of the data does not require higher- 
order modelling of the limb darkening). We use the Ohta et al. 
(2005) analytical description of the RM effect, and a Keplerian 
radial-velocity orbit. We compute the 3-D position of the planet 
relative to the star at each date. Two additional parameters are 



introduced to describe the photometric signal of the secondary 
eclipse: a planet-to-star surface-brightness ratio in the Spitzer 
band, and an out-of-transit continuum flux. The brightness vari- 
ation of the planet around the time of secondary eclipse is ne- 
glected. 

3.4. MCMC 

We use a chain with 200000 steps, starting with the parameters 
in M09. The size of the steps is set to 0.05 times the uncertain- 
ties quoted by that study. The free parameters are the follow- 
ing: P, Tq, i, e, m, V ( 5" 4 , M p \, R p i, M s ,R s ,dT, A, V rot , respectively 
the period, epoch of periastron, orbital inclination, orbital ec- 
centricity, argument of periastron, center-of-mass radial veloc- 
ity (for each spectrograph independently), planet mass and ra- 
dius, star mass and radius, planet-to-star surface brightness ra- 
tio in the Spitzer band, spin-orbit angle and projected stellar 
spin. We use Gibbs sampling for the steps in the chain (next 
step accepted with a probability equal to the ratio of likeli- 
hoods). We include the constraints from stellar evolution mod- 
els directly in the prior and in the merit function used in the 
MCMC: the MCMC chain moves in an interpolation of the 
Girardi et al. (2002) stellar evolution models with uniform steps 
in mass, age and metallicity (corresponding to a flat prior in 
these quantities). Each evolution model produces values for the 
stellar mass, radius and temperature, M s , R s and T e g, that are 
used to estimate the merit function. How this method solves 
Bayes' theorem is studied in more details in the context of stel- 
lar ages in Pont & Eyer (2004). 

3.5. Merit functions 

The merit function in the MCMC is taken as the likelihood 
of the observations given the model, assuming gaussian error 
distributions and purely "white+red" noise (Pont et al. 2006), 
using all the photometric and radial velocity data, as well as the 
spectroscopically determined stellar parameters: 



I 

phot 



(-^obs ^mod) 

°- 2 Fw +Ncr 2 Fr 



■Z 

rv 



(Vobs - V mo d) 
0-2 , at 2 

(S obs — S mod) 2 



(1) 

(2) 



where F are the photometric data, V the radial velocity 
data, S the stellar parameters, and the obs and mod subscripts 
denote observed and model values respectively. The w and r 
subscript indicate white and red noise parameters, and N is the 
number of data points during the correlation timescale. F mo d 
is the model light curve for the transit and secondary eclipse, 
V mo d the model radial velocity curve including the RM anomaly 
during transit, and S mod the set of observable parameters of 
the host star according to the Girardi et al. (2002) models. 
The sums are made respectively on the individual photometric 
observations, spectroscopic observations, and input stellar pa- 
rameters of the evolution models (age, mass and metallicity). 
We take the duration of the egress as the relevant correlation 
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timescale to evaluate N (~ 25 points for photometric time se- 
ries, 5 points for the radial velocity sequence). 

3.6. Prior distributions 

We define the steps in the MCMC corresponding to flat 
prior distributions in the following quantities: P, 7q, cos i, 
Vi, ecosm, esintzr, K^Jl - e 2 , R p /R s , M s , r s , [Fe/H], dT, 
cos A. The combinations were chosen either for physical rea- 
son (isotropic orientation of spin and orbit in space, random 
epoch), to avoid strong covariance in the MCMC (radius ratio 
instead of radius, eccentricity dependence of K integrated in 
the prior) or linear when the prior is not significant compared 
to the observational constraints (P, Vi). 

The prior distribution in the parameters of the parent star 
(M S ,R S , T e g) corresponds to a fiat distribution in age, mass and 
metallicity, as described above. We therefore do not reduce the 
information from stellar evolution models to a single term in 
the merit function, but make use of the full information pro- 
vided by the stellar evolution models, thus taking into account 
the correlation of stellar mass, radius and temperature, and the 
respective probability of different models for a field star in the 
solar-neighbourhood. 

We set the stellar rotational velocity to V ro , = 1.8 kms -1 
(Fischer & Valenti 2005). This parameter influences the shape 
of the spectroscopic transit radial velocity curve, but is poorly 
constrained by the current data. The results do not vary signif- 
icantly if we repeat the procedure with a moderate uncertainty 
on this value (up to about 1 kms _I ). If this parameter is left 
completely free, the Markov chain does not converge, because 
different combinations of the rotation velocity, spin-orbit an- 
gle and impact parameter produce similar predictions for the 
observed portion of the spectroscopic transit. 

Altogether, three prior distributions have a significant ef- 
fect on the solution: the orbital inclination angle (penalizing 
grazing transits compared to a non-Bayesian fit), the stellar age 
prior (penalizing rare or inexistant combinations of stellar pa- 
rameters), and the stellar rotation velocity (lifting the main de- 
generacy in the spectroscopic transit). 
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Fig. 4. Geometry of the HD 80606 system according to our best-fit 
solution, (a) from above the orbit, (b) seen from Earth. 
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Fig. 5. Posterior probability distribution for the projected spin-orbit 
angle. 



4. Results and discussion 

Table 2 gives the median values and central 68% confidence 
intervals for the parameters of the HD80606 system given by 
the MCMC integration. Note that because the median values 
correspond to different individual solutions for each parameter, 
the orbital parameters in the table do not correspond to the best- 
fit orbit. The best-fit orbit is P = 1 1 1.43605 days, K = 476.48 
ms" 1 , e=0.9332, T = 2454424.8529 BJD, w = 300.71°. 

Figures 5 and 6 show the probability distribution functions 
for the spin-orbit angle and transit duration (first to fourth con- 
tact). 

The strongest covariance between the output parameters is 
that of the impact parameter with the transit duration. In gen- 
eral, the transit duration depends mainly on the radius of the 
host star and impact parameter. When the whole transit is mea- 
sured, the main covariance is between impact parameter and 




15 
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Fig. 6. Posterior probability distribution for the transit duration (first 
to fourth contact) 
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Center-of-mass velocity V (Elodie) 
Vo(SOPHIE) 
AVo(Keck) 
AVo(HET) 

Orbital period P 
Orbital eccentricity e 
Velocity semi-amplitude K 
Epoch of periastron T 
Argument of periastron m 
Orbital inclination i 
Semi-major axis a 

Epoch of transit T u 
Epoch of eclipse T e 
Transit duration 7V 4 
Transit duration T 2 _ 3 
Radius ratio R p /R s 
Spin-orbit alignment A 
Impact parameter b 

Star Mass M s 
Star Radius R s 

Planet mass M p 
Planet radius R„ 



3.787 ±0.004 km s- 1 
3.910 ±0.004 km s" 1 
-0.001 ±0.002 km s" 1 
-0.019 ±0.003 km s-' 
11 1.4357 ±0.0008 days 
0.9332 ± 0.0008 
474 ± 4 ms" 1 
2 454 424.852 ± 0.008 BJD 
300.80 ± 0.22 
89.32 ± 0.06 
0.449 ± 0.006 AU 

2454876.316 ±0.023 BJD 
2454424.719 ± 0.009 BJD 
11.9 ± 1.3 hours 
9.6^ g hours 
0.103 ±0.003 
50° [14- 111]° 
0.75 ± 0.06 

0.97 ± 0.04 M 
0.978 ± 0.015 Ro 

3.94 ±0.11 Mj 
0.98 ± 0.03 Rj 



Table 2. Parameters for the HD 80606 system. Uncertainties from the 
68% central probability interval of the posterior distribution traced by 
the Markov chain. 



host star radius. However, in our case, the star radius is rela- 
tively well constrained by the duration of the secondary eclipse, 
and the transit duration is weakly constrained because the tran- 
sit ingress was not observed, so that the main covariance is be- 
tween the transit duration and impact parameter (higher impact 
parameters corresponding to shorter transit). Fig. 7 shows the 
posterior probability density in the D vs. b plane, and in the D 
vs R p plane. The core of the probability distribution in each pa- 
rameter is well described by Normal distributions, with a tail of 
low-probability solutions at high impact parameter, low transit 
duration and larger planet. 

Figure 8 shows the probability density in stellar mass and 
radius. The main correlation is M ~ R 3 , because the transit 
parameters are degenerate in these quantities. The sharper limit 
towards larger masses is due to stellar evolution models. 

4.1. Uncertainty intervals 

F09 quote very small uncertainties for the transit duration, or- 
bital inclinations and planetary radius (1.029 ± 0.017 Rj, 12.1 
± 0.4 hours, 89.285 ± 0.023 degrees). The uncertainties are 
small because F09 keep several important parameters fixed (the 
host star radius and orbital parameters), and purely uncorre- 
lated noise is assumed in the photometry. These simplifications 
are useful for a rapid initial analysis, but will yield underesti- 
mated uncertainties. As extensively discussed in the context of 
previous observations of transiting planets with high-cadence 
photometry, the uncertainty on the density (M/?~ 1/3 ) of the pri- 
mary and the systematics in the photometric noise are actually 



0.95 
Stellar Mass [A}.,; 



Fig. 8. Posterior probability density: stellar mass and radius. 



the dominant sources of uncertainties on most system parame- 
ters, including the planetary radius, inclination angle and transit 
duration. This explains why our uncertainty intervals are larger 
than those quoted in F09, in spite of being based on a larger 
ensemble of data. 

The larger uncertainties are also due to the fact that the 
M09, F09 and G09 light curves indicate slightly different 
shapes for the transit egress (M09 and G09 favour a longer 
egress, therefore a higher-latitude transit across a larger star). 
The shape difference between the three lightcurves clearly il- 
lustrates that correlated noise dominates the error budget. The 
presence of correlated noise is also apparent in the SOPHIE 
sequence, since no parameters for the RM effect can repro- 
duce the sudden jump shortly after the beginning of the data 
sequence. 

For these reasons, as well as the arguments and examples 
presented in Pont et al. (2006), the probability distributions de- 
scribed by the MCMC integration including correlated noise, 
uncertainties on all parameters, and host star properties con- 
strained by stellar evolution models, arguably provide a better 
description of the actual implications of the data in terms of 
physical parameters and confidence intervals. 



4.2. Spin-orbit angle 

The posterior probability distribution for the projected spin- 
orbit angle A corresponds to configurations in which the planet 
crosses mainly the receding limb of the star. An aligned orbit 
(A = 0) is excluded to a > 95% level of confidence. The MCMC 
calculation shows that a spin-orbit alignment in the system is 
made unlikely by the combination of the shape of the SOPHIE 
radial velocity time series during transit and the fact that the 
MEarth photometry excludes a low-latitude transit by placing 
a strict upper limit on the transit duration. 

The result is sensitive to the level of correlated noise as- 
sumed in the radial velocity data. This is easily understood: 
since the amplitude of the RM anomaly is of the order of 15 
ms -1 , assuming systematics of similar amplitude can recon- 
cile the observed shape with that expected in case of spin-orbit 
alignment by attributing most of the observed variation to sys- 
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Fig. 7. Posterior probability density: transit duration vs impact parameter (left) and transit duration vs planetary radius (right) 



tematics. The presence of some systematics in the SOPHIE 
data is apparent from the mismatch between the first few mea- 
surements during transit, to the level of a few ms -1 , but the 
excellent correspondence between the moment of photometric 
and spectroscopic egress indicates that the systematics do not 
dominate the SOPHIE sequence, and therefore that such a high 
value of <r r is unlikely (note that in our standard solution we 
already increase the correlated noise for the first measurements 
in the night to 12 ms _1 , which is conservative). 

Therefore, our conclusion is that the combined data in- 
dicate a tilted system, with the planet crossing the receding 
side of the star during the transit. HD 80606 is the second 
highly tilted planetary orbit known, after XO-3 (Hebrard et al. 
2008; Winn et al. 2009), out of the twelve measured systems 
(Fabrycky & Winn 2009). More strikingly, four of these plan- 
ets have eccentric orbits, and two of those have tilted orbit. We 
may therefore be seeing early indications that in addition to 
having a broad eccentricity distribution, extrasolar gas giants 
have a bimodal distribution of spin-orbit angles (Fabrycky & 
Winn 2009). A possible caveat is that both XO-3 and HD 80606 
are high-mass gas giants, which may have formed differently 
than Jupiter-mass planets (Ribas & Miralda-Escude 2007). 

4.3. Planetary radius 

The size of the planet is well constrained by the data, in spite 
of the fact that only the transit egress was observed. This may 
seem surprising. It is due to the observation of the secondary 
eclipse and the peculiarity of the orbit of HD80606b: the planet 
is much nearer to the star during the secondary eclipse than 
during the transit, because of the large eccentricity and the po- 
sition of the periastron. As a result, the mere fact that the planet 
is transiting implies that the secondary eclipse occurs at low 
latitude behind the star (see Fig. 4). This lifts the usual degen- 
eracy between stellar radius and impact parameter, and means 
that the former can be derived reliably from the duration of 
the secondary eclipse, even in the absence of the constraint 
from the fractional duration of the ingress/egress. The proba- 
bility distribution of the stellar radius is further narrowed by 
the Bayesian prior on the inclination angle, which penalizes 



high-latitude transits, by the shape of the transit egress, and by 
stellar evolution models. 

The main remaining uncertainty on the size of the planet is 
due to global systematics in the photometry. The slightly longer 
egress duration favoured by the M09 and G09 data would indi- 
cate a higher impact parameter, therefore a larger star and a big- 
ger planet. However, this is partly compensated by the fact that 
these datasets also favour a slightly shallower transit, therefore 
a smaller planet. Altogether, the uncertainty on the planetary 
radius is small enough for useful comparison with models. 

We confirm that HD 80606 is not a "bloated" hot Jupiter, as 
discussed in M09, and has a size compatible with current mod- 
els for irradiated gas giants. Given its large mass, a planet like 
HD 80606b would require more additional energy that Jupiter- 
mass planets to be inflated to a higher radius. However, two 
other exoplanets with similar masses, CoRoT-Exo-2b (3.3My) 
and OGLE2-TR-L9b (4.5M y ) have radii around 1.5 Rj, show- 
ing that such inflated sizes are indeed possible for these heavy 
planets. The two other known transiting planets in the 3-5 Mj 
mass range, HD 17156b and WASP-lOb, have radii of 1.0 and 
l.lRj respectively. 

The size of these planets is conspicuously correlated with 
the amount of flux received from the parent star. CoRoT-Exo- 
2b and OGLE2-TR-L9b follow close orbits (1.7 and 2.4-day 
periods), while HD 17156b and HD 80606b have the widest 
orbits of known transiting planets, with WASP- 10b being an 
intermediate case. The observed correlation between size and 
incident flux for transiting gas giants may therefore extend to 
several Jupiter masses, which provides another benchmark that 
any successful explanation of anomalous exoplanet radius must 
meet. 

Figure 9 shows the evolution of the radius of HD 80606b 
according to models described in Guillot et al. (2006), with 
the position of our reference solution indicated. The differ- 
ent panels show the radius evolution of a planet with a solar- 
composition envelope (containing about 25 M ffi in heavy ele- 
ments), but either without core (left panel), with a 100M e core 
(middle panel) or with a 200M ffi core (right panel). Clearly, al- 
though uncertainties in the distribution of heavy elements and 
equations of state have to be taken into account (e.g. Baraffe 
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et al. 2008), our models point towards the existence of a large 
mass in heavy elements (at least 60 M ffi ) present in the planet. 
All possibilities above 60 and 200 e are consistent with the mea- 
sured radius. Solutions with even larger cores are possible, but 
we regard them as unlikely given that realistic critical core 
masses are smaller than 100M ffl (Ikoma et al. 2006) and plan- 
ets with masses larger than Jupiter tend to scatter planetesimals 
much more efficiently than they accrete them (e.g. Guillot & 
Gladman 2000). 

Given that HD80606 is metal-rich, this large inferred core 
mass is consistent with a correlation between heavy elements 
in the star and in the planet (Guillot et al. 2006; Burrows et al. 
2007; Guillot 2008). It should be noted that if the planet had the 
same composition as the star, its total mass of heavy elements 
should be of order 50 M e , clearly smaller than the estimates 
that are derived here. 

Interestingly, using the concept of a maximum mass in 
heavy elements in the planet yields an upper limit on the rate 
of tidal dissipation that it may undergo. As shown by the dotted 
lines in Fig. 9, a maximum rate of dissipation compatible with 
the observations and models with core of at most 200 M ffl is of 
order 10 27 erg s _1 . While this may seem large, it is several times 
smaller that the dissipation due to tides for this planet if it has 
been brought in by a Kozai mechanism (Wu & Murray 2003), 
estimated by GM t M p /(a p T migmtioa ) w 5 x 10 27 ergs~' (where 
we used a p as the present semi-major axis of the planet, and 
^migration ~ 1 Ga from WM03). Several possibilities exist: (i) 
The planet may have a tidal Q larger than the 3 x 10 s assumed 
by WM03 so that it would have migrated on a longer timescale; 
(ii) Heat dissipation may occur in a shallow region close to the 
atmosphere of the planet and be reradiated quickly after the ap- 
proach to the star; (iii) The planet may have a core that is much 
larger than envisioned in the present study. Except in the last 
case, our calculations indicate that the planet's intrinsic effec- 
tive temperature should be of order 200 to 300 K, much smaller 
than the 600 to 700 K envisioned by Wu & Murray (2003) and 
Laughlin et al. (2009). 



4.4. Future observations 

Given the brightness of the host star and its peculiar orbit, HD 
80606 is likely to remain an important target for further studies. 
Table 3 shows the predicted times of ingress and egress for the 
next transit events according to our reference solution and 68% 
central confidence intervals. 



Date 


ingress [BJD] 




egress [BJD] 




5 June 2009 


2454987.53 ± 


0.05 


2454987.974 


± 0.005 


24 September 2009 


2455098.96 ± 


0.05 


2455099.409 


± 0.005 


13-14 January 2010 


2455210.40 ± 


0.05 


2455210.844 


± 0.006 


5 May 2010 


2455321.84 ± 


0.05 


2455322.281 


± 0.006 


24-25 August 2010 


2455433.27 ± 


0.05 


2455433.717 


± 0.006 


14 December 2010 


2455544.71 ± 


0.05 


2455545.152 


± 0.007 


4-5 April 2011 


2455656.14 ± 


0.05 


2455656.588 


± 0.008 
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Table 3. Predicted ephemerides for the next transit events 



2468 2468 2468 
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Fig. 9. Evolution of the planetary radius according to models for HD 
80606b. Left: without heavy-element core. Middle: with a 100 M e 
core. Right: with a 200 M e core. From bottom to top, the lines corre- 
spond to different assumption on the physics of the evolution models 
used to match the sizes of transiting planets: blue: standard models 
with no extra heating; orange: dissipating 0.5% of incoming stellar 
heat flux energy in the interior of the plane ; red: opacities increased 
by a factor 30. The dotted (blue) lines correspond to models with an 
extra luminosity dissipated at the center of the planet L tide = 10 26 and 
10 27 erg s , respectively. 



When the transit duration will be measured by observation 
of future transit events, the corresponding value of the most 
probable impact parameter and planetary radius can be read 
directly off Fig. 7. 

5. Conclusion 

Our analysis of combined photometric and spectroscopic ob- 
servations of HD 80606 have led to new estimates of the stel- 
lar and planetary parameters of the system. These values allow 
relatively tight constraints on the planet's composition, with a 
likely mass in heavy elements between 60 and 200 M®, which 
imply a rather efficient accretion of planetesimals most proba- 
bly during the planet's formation and early evolution. 

We have also shown that the planet's orbit is probably not 
aligned with the star's spin: this strengthens the argument that 
the planet may have migrated inwards through a Kozai mech- 
anism and thus acquired its large eccentricity (Wu & Murray 
2003). 

Our analysis can be readily put to the test with future ob- 
servations of planetary transits, as we predict a total duration 
of the eclipse of T1-4 = 11.9 ± 1.3 hours. Furthermore, the 
analysis of the planet's primary and secondary transits, and of 
the planet's irradiation in the infrared when far from the star 
will also be extremely fruitful to understand how tidal energy 
is dissipated: we predict that the intrinsic effective temperature 
of the planet is between 200 and 300K, smaller than the zero- 
albedo mean equilibrium temperature due to stellar irradiation, 
400K. In contrast, larger effective temperatures (~700K) have 
been hypothetized, due to efficient dissipation of heat in the 
planet (Laughlin et al. 2009). 
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